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Systems that exhibit pattern formation are typically driven and dissipative. In the early universe, 
parametric resonance can drive explosive particle production called preheating. The fields that are 
populated then decay quantum mechanically if their particles are unstable. Thus, during preheating, 
a driven-dissipative system exists. We have shown previously that pattern formation can occur in 
two dimensions in a self-coupled inflaton system undergoing parametric resonance. In this paper, we 
provide evidence of pattern formation for more realistic initial conditions in both two and three di- 
mensions. In the one-field case, we have the novel interpretation that these patterns can be thought 
of as a network of domain walls. We also show that the patterns are spatio-temporal, leading to a 
distinctive, but probably low-amplitude peak in the gravitational wave spectrum. In the context of 
a two-field model, we discuss putting power from resonance into patterns on cosmological scales, in 
particular to explain the observed excess power at 100/i _1 Mpc, but why this seems unlikely in the 
absence of a period of post-preheating inflation. Finally we note our model is similar to that of the 
decay of DCCs and therefore pattern formation may also occur at RHIC and LHC. 
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I. INTRODUCTION 

Much recent work has been devoted to the topic of pre- 
heating in inflationary cosmology. Preheating is a stage 
of explosive particle production which results from the 
resonant driving of particle modes by an inflaton oscil- 
lating in its potential at the end of inflation . 

In regions of parameter space where parametric res- 
onance is effective, much of the energy of the inflaton 
is transferred to bands of resonant wave modes. This 
energy transfer is non-thermal and can lead to interest- 
ing non-equilibrium behavior. Two examples of the non- 
equilibrium effects that can be produced are non-thermal 
phase transitions and baryogenesis The non- 

thermal phase transitions induced during preheating can 
sometimes lead to topological defect formation [p~2|— p~5| , 
even at energies above the eventual final thermal tem- 
perature. Furthermore, non-linear evolution of the field 
when quantum decay of the resonantly produced parti- 
cles is negligible leads to a chaotic power-law spectrum 



of density fluctuations |16,fL7|. 

In a previous letter [19|, we presented evidence for a 
new phenomenon that can arise from preheating: pattern 
formation. It has long been known that many condensed 
matter systems exhibit pattern formation^. Examples of 
pattern forming systems which have been studied are rip- 
ples on sand dunes, cloud streets and a variety of other 
convective systems, chemical reaction-diffusion systems, 
stellar atmospheres and vibrated granular materials. All 
of these physical systems have two features in common. 
They are all driven in some manner, i.e. energy is input 
to the system, and they are all dissipative, usually be- 
ing governed by diffusive equations of motion. Typically, 
patterns are formed in these systems in the weakly non- 
linear regime before the energy introduced into the sys- 
tem overwhelms the dissipative mechanism. Sometimes, 
patterns persist beyond the weakly non-linear regime as 
well. 

At the end of inflation, the inflaton <f> is homogeneous 
and, in most commonly studied models, oscillating about 
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'For an extensive review of pattern formation in condensed 
matter systems, see |2l[j22|]. 
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the minimum of its potential. This oscillation gives an 
effectively time-dependent mass to fields which are cou- 
pled to the inflaton, including fluctuations of the infla- 
ton itself. The time-dependent mass drives exponen- 
tial growth in particle number in certain bands of wave 
modes. However the fields into which the inflaton can de- 
cay resonantly are also unstable to quantum decay. For 
these reasons, at the end of inflation, we are considering 
fields which are driven, due to resonant particle creation, 
and also dissipative, due to quantum decay. In |l9|], we 
were able to show that pattern formation occurs in a 
chaotic inflationary model with a self-coupled inflaton in 
the weakly non-linear regime. 

In this paper, we primarily consider a A0 4 theory with 
the addition of a phenomenological decay term to mimic 
the inflaton's quantum decay. This model without the 
decay term has been studied extensively in the literature 
|,|jl|Jl|Jl^||,||, and a similar model including the 
decay term has also been studied |10|. 

In 1 19 1, we used restricted initial conditions. We only 
seeded the resonant band with small fluctuations of or- 
der ~ 10~ 3 of (4>), then simulated the field's evolution. 
We found that the resonant modes interacted and formed 
patterns. Here we use initial conditions appropriate to 
the vacuum at the end of chaotic inflation. We also ex- 
tend our study to a 3-dimensional volume. In fl9|] , we did 
not point out the spatio-temporal nature of the patterns, 
which was not evident at the time due to an unfortunate 
coincidence in the form of the pattern at the timesteps 
at which we viewed the data. Here we note the spatio- 
temporal behavior. We also discuss resonance giving rise 
to patterns on cosmological scales. 

The paper is ordered as follows: In section O, we 
present the model we are investigating. In section p|, 
we discuss the initial conditions appropriate for the end 
of inflation. In sections IV, we present the results of our 
simulations in two and three spatial dimensions. In sec- 
tion 0, we extend our analysis to the two-field case and 
discuss the possible implications of our results. This is 
followed by the conclusions. 



II. THE X(j) A MODEL WITH 
PHENOMENOLOGICAL DAMPING 

Our field equation in comoving coordinates is 



3H <j> + 70 - 



1 



+ A0 3 = 



(1) 



where 7 is a decay constant, A is the self-coupling of the 
field and H = a/ a. We convert to conformal time dt — 
a{r)dr and introduce the field (p = a<j). Upon rescaling: 
t — > r/^/Xipfi, x — > x/V\ipR and tp — > pp>R, where 
subscript R denotes "at the start of reheating" , we obtain 
a new equation 



ip" + aTp - VV- 0'r+ — )<p + p 3 = 0, (2) 

a 

where T = j/y/XipR. Further simplification is possible if 
we note that in A</> 4 theory, the averaged equation of state 
during preheating is that of radiation, therefore a" ~ 0. 
We use a = 1 + Hrt. 

It should be noted that pattern formation in the infla- 
ton system is conceptually distinct from condensed mat- 
ter systems for at least two reasons. First, the equations 
we study are wave equations with damping, not diffu- 
sive equations. Secondly, we expect wave patterns to be 
formed while the homogeneous mode decays, therefore 
pattern formation will be a temporary phenomenon, at 
least in the model above in which gravity is neglected. 
The driving in Xcj) 4 preheating comes from the large ini- 
tial value of the inflaton at the end of inflation, causing 
the field to roll and oscillate in its potential. This should 
be considered in contrast to the typical condensed mat- 
ter system, in which energy is introduced via boundary 
conditions (in a convective system) or by a vibrating bed 
(in a granular material system), and the energy input is 
essentially constant. 

At the end of inflation, we may expand tp about 
a homogeneous piece and then linearize Eq. (Q). Let 
ip(r, x) = $(r) + ip(r,x), where j ipd 3 x — 0. We ob- 
tain 



$" + $3 = Q 

K + aI >k + ( fc2 - oT + 3$ 2 (t)) Vk = 0, 



(3) 
(4) 



where we have taken the Fourier transform in the latter 
equation. The solution to Eq. (0) is $(r) = cn(r; l/y/2), 
and therefore, for T = 0, Eq. (|4|) becomes a Lame equa- 
tion. 

For r = 0, the resonant modes lie in the interval |23l 



- < e < vs, 



(5) 



and have an amplitude of the form exp(/Xfer), where \Xk 
is the characteristic exponent or Floquet index. It should 
be noted that the wavelengths of the resonant modes are 
of order the Hubble radius immediately after inflation. 

For r ^ 0, we can introduce ip = tp exp(i J aTdr) 
giving 



k 2 



-a'T 
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a *r' + 3&{T)]ifr t = 0. (6) 



Now the resonance band becomes time-varying. The 
modes which can be in resonance at any moment in time 
satisfy | < k 2 ^ < \/3, where 



^eff 



Vr 



: a 2 T 2 . 



(J) 



However the modes we are actually interested in are 
■0k/a, and these will only grow for times r < 4\nk\/ HrT . 
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This is because we typically have r <C < Hr ^ 1 (in 
rescaled units). Thus, even in the absence of backreac- 
tion||| a combination of the expansion of the universe and 
quantum mechanical decay of the inflaton serves to take 
the system out of resonance. During this time the values 
of k 2 in resonance go from being shifted upwards by an 
amount T to an amount Effectively this means the 
resonance band is smeared out in /c-space. 

While the above perturbative analysis is helpful it is 
suitable only for early times as resonance will soon take us 
away from the linear regime. The effect that we are trying 
to isolate is intrinsically non-linear so we now resort to 
numerical simulation of the field equation (En. 



III. INITIAL CONDITIONS 

The initial conditions for $ are those of the end of slow- 
roll, which is normally supposed to be when a = 0. It is 
sufficient for our purposes to set $(0) = 1 and $'(0) = 0. 

The fluctuations in the inflaton are quantum in origin. 
Super-Hubble modes exist which were generated during 
inflation and have remained fixed in amplitude after leav- 
ing the Hubble volume. However these modes do not un- 
dergo amplification in this model, so it is not important 
to quantify them precisely. Of more significance for us 
are the sub-Hubble modes at the end of inflation. At 
this time V is in fact zero and it is appropriate to con- 
sider any field to be in the vacuum state. Then the initial 
conditions for the tp^ are given by the usual results: 



|A k |e 
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(8) 



where uj\ — k 2 + 3<f>(0) 2 , is a number randomly taken 
from a Gaussian distribution with zero mean and unit 
variance, and is a random number taken from the 
interval [0, 1]. 



IV. SIMULATION RESULTS 

A. Two dimensions 

To simulate the evolution of the field, we discretize the 
spatial derivatives to fourth-order in Ax, and we use a 
leapfrog integrator which is accurate to second order in 
At. 

Setting the box size to 256 gridpoints per dimension 
and such that the resonant wave number in the box is 



16 is enough to give many different resonant modes in 
the box, but still have good resolution of the wave, so 
this is the box size we used. We use periodic boundary 
conditions. 

We set out to identify the weakly non-linear regime. 



Khlebnikov and Tkachev |16| showed that, without the 
decay terms, the self-coupled inflaton system's non-linear 
time evolution proceeds as follows: First, the resonant 
band amplitude grows. Next, when the amplitude in 
the resonant band is high enough for non-linear effects 
to become important, period doubling occurs and sub- 
sidiary peaks develop in the power spectrum. Further 
peaks then develop and the spectrum broadens and ap- 
proaches an exponential spectrum. We tuned T such that 
the amplitude of the resonant band grew, but little period 
doubling occurred. In this regime, only resonant mode 
wavelengths exist in the box and they interact with each 
other non-linearly. 

In the expanding case, the increase of the effective 
damping coefficient r e g- = aT with increasing scale fac- 
tor a makes it easier to keep the system in the weakly 
non-linear regime compared to the non-expanding case. 

We found the smallest value for T such that the system 
remained weakly non-linear was of order 10 -5 . This is 
two orders of magnitude smaller than that in the non- 
expanding case. 

In Figure [l] we plot a superposition of the power spec- 
trum at various times during a simulation with T = 
5 x 1CP 5 . It is possible to see that the system stays in 
the weakly non-linear regime for the entire simulation. 



0.012 • 
0.01- 



25 50 75 100 125 150 175 200 
FIG. 1. A superposition of power spectra taken from a sim- 
ulation with r = 5 x 10 . We plot the amplitude of the 
power spectra vs. wave number. Note that the period dou- 
bling modes are only weakly populated, indicating that the 
simulation is in the weakly non-linear regime. 



^It should be noted that backreaction profoundly alters the 
picture that emerges from a perturbative treatment; see e.g. 
[ p9[ . In our case backreaction is necessary but also necessarily 
small. 



When wave patterns form, the specific pattern which 
arises is due to the non-linear interaction of the wave 
modes. The amplitude of wave modes separated by dif- 
ferent angles grows at different rates. Modes separated 
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by angles with the fastest growing amplitudes dominate 
the solution and form the wave pattern. 

It is also possible that patterns vary temporally. This 
is the situation we find in the A^> 4 model. What we 
see in both the expanding and non-expanding cases is 
that a pattern emerges from the fluctuation background, 
then the peaks and valleys begin to move relative to each 
other. 



Initial conditions for a two-dimensional 

5 



FIG. 2. t = 

simulation with F = 5 x 10 -5 . The magnitude of <j>(k x , k v ), the 
Fourier transform of the inflaton, is indicated by the shading: 
brighter regions correspond to regions of larger amplitude. 
The zero mode has been deleted for plotting purposes, and 
the surrounding modes populated with vacuum amplitudes 
from Eq. H. Only the region of interest is plotted. 




FIG. 4. t = 1020. The final wave pattern in Fourier space. 

The temporal dependence is almost periodic. Peaks 
and troughs in the field energy align along one direction 
in a ripple-like pattern, then the pattern flips to align in 
a direction orthogonal to the original direction. We say 
the dependence is 'almost' periodic because the field flips 
back and forth, but the timing of the flips varies as the 
field evolves. This leads us to believe that, for instance, if 
there were a background driving field that gave constant 
energy input (as opposed to the decaying background in 
chaotic inflation) the flipping would be truly periodic. 




FIG. 5. t = 0. The initial conditions in configuration space. 
The inflaton 4>(x,y) is plotted, again with brighter regions 
corresponding to larger values of </>. 



FIG. 3. t — 840. After resonance begins to boost the ampli- 
tude of the resonant mode. Notice the brightening (increasing 
amplitude) of the resonant modes. 
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a'T) 



(9) 



FIG. 6. t = 840. Wave pattern at intermediate time. 



in a radiation dominated universe. The barrier height is 
(aT/2) 2 which is constant in time. Thus dissipation not 
only leads to an overall damping of the field, but also 
leads to degeneracy in the minima of the effective poten- 
tial. This suggests a novel interpretation of pattern for- 
mation in this model: the pattern is actually a network of 
domain walls separated by the characteristic wavelength 
of the resonance. If T increases then the barrier height in- 
creases, reflecting the fact that, the non-linear coupling of 
the modes necessary for pattern formation is suppressed. 
On the other hand, for small T, the barrier is not high 
enough to prevent the field from probing both minima, 
and again patterns do not form. This is because the sys- 
tem has become strongly non-linear. 



B. Three dimensions 




FIG. 7. t = 1020. Wave pattern at the end of the simula- 
tion. 

In Figures |[ || and ^ we present snapshots of the evo- 
lution of the Fourier transform of the inflaton in two di- 
mensions in an expanding universe. And in plots § § 
and [7] we present snapshots of the evolution of (f)(x, y) 
in configuration space in two dimensions in an expand- 
ing universe. Notice the change in direction between the 
ripples in Figure ^ and |^. 

In the expanding case (in contrast to the non- 
expanding case) , the patterns are more clearly delineated 
and look less noisy. This is due to the fact that the field 
equation has an effective symmetry breaking potential 
and this leads to a restoring force on the fluctuations in 
ip. From Eq. (0) we deduce 



Using the same simulation techniques described above, 
we also find spatio-temporal pattern formation in three 
dimensions. We plot the field in configuration space in 
Figures |[ || and 10, again with r = 5 x 10~ 5 . The behav- 
ior of these patterns is similar to those found in two di- 
mensions: the pattern forms at the resonant wavelength, 
then the peaks and valleys begin to move with respect to 
each other. In these simulations, we set the physical box 
size such that there were eight resonant wavelengths per 
dimension in the box. Overall there were 64 gridpoints 
per dimension and we neglected expansion in these sim- 
ulations. 




FIG. 8. t = 0. Initial conditions for a three-dimensional 
simulation with T = 5 x 10 -5 . Here 4>(x,y,z) is plotted, 
again with brighter regions indicating larger field values. 
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FIG. 9. t = 100. Wave pattern at intermediate time. 




FIG. 10. t — 300. Wave pattern at the end of the simula- 
tion. 



V. DISCUSSION 

A. Gravitational waves 

Since the patterns that we see vary in time, we expect 
there to be a peak in the gravitational wave spectrum at 
the resonant frequency for preheating in the weakly non- 
linear regime. Similar peaks in the gravitational wave 

in simulations of un- 



spectrum have been seen 24 21 
damped preheating. These peaks correspond to the res- 
onant and period doubled frequencies. In the pattern 
forming regime, we expect to see similar peaks of lower 



amplitude, since the dissipative terms keep the resonant 
peak in our simulation at lower amplitude than the un- 
damped system. The peak in the gravitational wave spec- 
trum from pattern formation will be of smaller amplitude 
than that of the fully chaotic system but since the pat- 
tern has directionality, this might help in extracting a 
signal. 



B. Large scale structure formation 

We would also like to discuss the possibility of patterns 
occurring on cosmological scales, since, if we could put 
resonant power at 100ft. _1 Mpc, we would have an expla- 
nation for the observed excess power found at these scales 
|p6| . Although such a scale is outside the Hubble ra- 
dius at the end of inflation, recently a number of authors 
p?| p3| have investigated the question of whether reso- 
nance can amplify modes with wavelength larger than 
the Hubble radius. It turns out that this is possible be- 
cause of the large-scale (many Hubble volumes) coher- 
ence of the inflaton at the end of inflation; super-Hubblc 
mode amplification can be thought of as down-scattering 
from the oscillating zero mode. However in the usual pa- 
rameter region of some models J30|,[3l],[}3| , this amplifica- 
tion does not seriously affect the post-inflationary power 
spectrum. This is due to the fact that during inflation 
in these models, fluctuations in matter fields are sup- 
pressed by a factor of a -3 / 2 compared to fluctuations 
in the inflaton. Therefore, during preheating, although 
super-Hubble modes are amplified, they cannot be am- 
plified enough to be significant. 

However consider the following model for preheating: 



(10) 



For q = g /A ~ 1, fluctuations in x are n °t suppressed 
during inflation |32|. This is because the effective mass 
of x is g<fi, and this is much less than H during slow roll. 

If we now suppose \ decays quantum mechanically into 
other particles, then we are led to a phenomenological 
equation of motion akin to Eq. ([!]) 



X + 3Hx + 1 X X 



1 



V 2 



x ■ 



2 ±2 

1 <t> X 



0. 



(11) 



With the same simplifications used earlier, and intro- 
ducing X — xaexp(i J draT x ), we obtain the the fol- 



lowing equation for the Fourier modes of X 



off 



+ <7$ 2 (r))X k = 0, 



(12) 



where now fc 2 ff = k 2 - 3aT x /2 - a 2 r 2 /4. This differs 
essentially from Eq. (^) because of the appearance of q. 
A partial plot of the instability bands in the (k 2 ^, q)- 
plane is shown in Figure [0]. When r x ^ we have 
a new feature: for modes with k small enough, we can 
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have fc 2 g- < 0, i.e. the system can probe the lower half of 
instability plot. 




FIG. 11. A portion of the instability plot for q ~ 1 in 
Eq. ([l2]). Darker shading indicates larger values of the char- 
acteristic exponent /ifc. The narrow resonance band of Eq. ([?]) 
is not indicated. 

If fcgk y is the physical scale today on which we want 
to have a feature in the power spectrum, then this corre- 
sponds to a co-moving mode kn at reheating such that 



7~Lr 



H an Hr ' 



(13) 



where Ti = a' /a = aH. For a bump at 100/i _1 Mpc, we 
find we need k^/Ti-R ~ 10~ 24 . In terms of our dimension- 
less quantities this becomes fc# ~ 1CP 24 which we want 
to be the principal wavelength in resonance. It follows 
that we have either k 2 ^ < or T x < 10 8 , but from our 
simulations the latter is ruled out because it will be too 
small to damp resonance and lead to pattern formation. 

If we consider Figure 11 it is readily apparent that it 
will be difficult to introduce a feature at the required 
scale. For g < 1 we can have resonance at k 2 = 10~ 24 
but also, ruinously, at all scales down to k = 0. In- 
deed the characteristic exponent becomes larger as k gets 
smalleif*| We can avoid this problem for q > 1 but only 
if q is fine tuned to an unacceptable degree. In any case 
the resonance band is broad here (this problem is further 
exacerbated in our situation by the fact that quantum 
mechanical decay tends to smear the resonance bands) 
and peaks for modes with k 2 q > 0. Thus it is unlikely 
that super-Hubblc mode amplification in this model can 
lead to a feature at 100/i -1 Mpc. 

However we now mention that it may not be necessary 
to require such amplification in order to have an effect on 



large scale structure. It is possible that a posi-preheating 
period of inflation |34|,[I],[35| occurs (or even several) and 
this may bring sub-Hubble modes, previously amplified 
in narrow band resonance, up to scales appropriate for 
structure formation. The growth in the scale factor will 
also be larger if we imagine the inflaton decays into a 
large number of fields Xii though we have not considered 
pattern formation in such models. 

It should be noted that T x is not constrained in the 
same way is. The latter has to be small so that radia- 
tive corrections do not spoil the required flatness of the 
inflaton potential. In the rescaled units used here, it is 
not hard to show Q that < 10~ 13 , where we have 
taken A ~ 10 -13 which leads to the correct magnitude 
of density perturbations from inflation. Thus our earlier 
results were obtained with a T which was too large to be 
physical. However the results should go through for the 
two-field model considered above. 



C. Patterns in heavy ion collisions 

As a final comment, we would like to point out that 
the system of equations which we have investigated is 
similar to that for heavy ion collisions. It is hoped that 
in coming experiments at RHIC and LHC energies will 
be obtained which are in excess of the critical temper- 
ature for the QCD chiral phase transition. The situa- 
tion may be modelled by an 0(4)-symmetric theory for 
the scalar fields <& = (a, ff) (37| , where one initially has 
(a) = (7?) =0. Due to the expansion of the plasma, the 
temperature quickly falls below the critical temperature, 
the 0(4)-symmetry is broken, and the system evolves to 
new equilibrium values (a) 0, (7?) = 0. It is during 
this time that "disoriented chiral condensates" (DCCs), 
domains in which the pion field develops a non-zero ex- 
pectation value in a certain direction, may be produced. 
The subsequent decay of a DCC gives rise to a charac- 
teristic signal in the detectors. 

The important point for us is that during DCC de- 
cay, long-wavelength pion modes are resonantly ampli- 
fied when a oscillates about the minimum of its effective 
potential |pq| . Furthermore a decays via a — > 2ir and 
it has been argued that this process can be modelled 
phenomenologically by the addition of a time-dependent 
friction term |39[| . Although some questions remain con- 
cerning the applicability of such a term, it is intriguing 
that a driven-dissipative system exists for DCC decay. 
This suggests that patterns might turn up at RHIC and 
LHC. We are currently working on the experimental sig- 
natures of such events. 



**It follows that, for 5 < 1, the quantum mechanical decay 
of the inflaton leads to a red tilt in the power spectrum. 
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VI. CONCLUSIONS 

In this paper, we have presented new evidence that 
there is a pattern forming regime in a A0 theory during 
preheating at the end of inflation. We have used vac- 
uum initial conditions appropriate to sub-Hubble modes 
at the end of inflation. We have shown that patterns 
arise in both two- and three-dimensions, and that, in 
the one-field case, they may be thought of as a network 
of domain walls. Furthermore since the patterns vary 
spatio-temporally, gravitational waves will be produced. 
However, relative to the gravitational waves produced in 
an undamped model, their amplitude will be small, and 
therefore, unless the directionality of the pattern can aid 
in detection, extremely difficult to detect. 

We have speculated on the possibility, in the context of 
a two-field model, of putting a resonant band at cosmo- 
logical scales, in particular at 100/i _1 Mpc. We conclude 
that this is not feasible because it is impossible to am- 
plify only the required modes. However this conclusion 
may change in a model which allows a subsequent period 
of inflation. 

Finally, we point out that heavy ion collisions at RHIC 
and LHC may be another physical system which exhibits 
pattern forming behavior. 
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